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Abstract 

We compare Monte Carlo results for the interface tension and interface energy 
of the 3-dimensional Ising model with Pade and inhomogeneous differential approx- 
imants of the low temperature series that was recently extended by Arisue to 17**^ 
order in m = exp(— 4/3). The series is expected to suffer from the roughening sin- 
gularity at M ~ 0.196. The comparison with the Monte Carlo data shows that the 
Fade and inhomogeneous differential approximants fail to improve the truncated se- 
ries result of the interface tension and the interface energy in the region around the 
roughening transition. The Monte Carlo data show that the specific heat displays a 
peak in the smooth phase. Neither the truncated series nor the Pade approximants 
find this peak. We also compare Monte Carlo data for the energy of the ASOS 
model with the corresponding low temperature series that we extended to order u^"^. 
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1 Introduction 



Properties of the interface separating coexisting phases of the 3-dimensional Ising 
model have found continuous interest in the hterature. 

The Ising interface undergoes a roughening transition at an inverse temperature 
jSr = 0.4074(3) [H that is nearly twice as large as the bulk transition coupling 
= 0.221652(3) 0. 

The roughening transition is believed to be of the Kosterlitz-Thouless (KT) 
nature 0, with the surface free energy having an essential singularity at f3r of the 
type 

/sing~exp[-A(/5-/5,)-i/2]. (1) 

Though the free energy and all its derivatives with respect to (3 stay finite at the 
roughening point, one has to expect that low temperature series for interface prop- 
erties suffer from the transition. 

The first low temperature expansion of the 3D Ising interface tension a was given 
by Weeks et al. to 9**^ order in the variable u = exp(— 4/3). Shaw and Fisher [Q ana- 
lyzed the series with the help of Pade and inhomogeneous differential approximants. 
They claimed that the Pade and differential approximants allow to compute the 
surface tension accurately for temperatures below the roughening point. 

Recently, Arisue put forward the series to 17**^ order in u |p. It is interesting to 
note that the coefficients of the series change their sign at order 13. This behavior 
does not come unexpected and confirms that the roughening transition of the Ising 
interface is of Kosterlitz-Thouless type: Expanding eq. (|I]) in the variable u, one 
also obtains a series with coefficients that change their sign at a certain order. The 
order where the change of sign happens depends on the non-universal parameters A 
and Pr- 

A completely independent method to study the Ising interface is the Monte 
Carlo method 0. For recent Monte Carlo work on the 3D Ising interface see, e.g., 
0, |, 1^, |l^, [Tl|, [12| and references cited therein. 

In we reported on a numerical study of properties of the Ising interface over 
the whole range from low temperatures up to the bulk critical point. In particular 
we determined the surface energy and, by integration over also the surface free 
energy and surface tension. 

In this paper, we give a more detailed account of the numerical results and com- 
pare them with Pade and inhomogeneous differential approximants for the extended 
low temperature series. In order to demonstrate that the disagreements of series 
and numerical data are not due to finite size effects, we provide data for various 
lattice extensions and demonstrate that the systematic errors in the determination 
of surface energy and surface tension are under control. 

The ASOS (absolute value solid-on-solid) model in two dimensions is an approx- 
imation of a phase separating surface in a 3-dimensional Ising model. The approx- 
imation neglects overhangs of the surface and fiuctuations of the bulk phases. For 
temperatures being low compared to the bulk critical temperature of the Ising model 
this approximation should describe the Ising surface quantitatively rather well. We 
extended the existing low temperature series for the ASOS model to 12*^ order and 



compared its Pade approximants with Monte Carlo results that we obtained with 
the VMR cluster algorithm [13|. 

This paper is organized as follows. In section ^ we introduce the notation for the 
3D Ising model and its ASOS approximation. We also define the interface tension 
and discuss the problem of how to obtain it from simulations of finite lattices. The 
algorithms used for the Monte Carlo simulations are specified in section |^. In section 
§ we quote the results and compare them with extrapolations of low temperature 
series. Conclusions are given in section 1^. 



2 The Models and Quantities Studied 
3-Dimensional Ising Model 

We consider a simple cubic lattice with extension L in x- and ?/-direction and with 
extension t in 2;-direction. The Ising model is defined by the partition function 

Z= ^M-PH). (2) 

The Hamiltonian if is a sum over nearest neighbor contributions, 

H = - J2 f^ij^i^j ■ (3) 

<i,j> 

The interaction energy is normalized such that (3 = 3/ k-gT, where k-g denotes Boltz- 
mann's constant, J is the interaction energy, and T is the temperature. 

The lattice becomes a torus by defining that the uppermost plane is regarded 
as the lower neighbor plane of the lower-most plane. An analogous identification is 
done for the other two lattice directions. For the Ising spins a we use two different 
boundary conditions: Periodic boundary conditions are defined by letting kij = 1 
for all links < i,j > in the lattice. To define antiperiodic boundary conditions 
in 2;-direction, we also set kij = 1 with the exception of the links connecting the 
uppermost plane with the lower-most plane. For these links we set kij = —1. 

Let us define 

G = -(lnZ„-lnZp) , (4) 

where the subscript a (p) means antiperiodic (periodic) boundary conditions. The 
surface tension can be defined as the limit 

G 

a = lim lim — . (5) 

With numerical simulations only a rather limited range of L and t values is accessible. 
Hence a careful discussion of finite size effects is needed. Let us express the partition 
functions of the periodic and antiperiodic Ising system in terms of the transfer matrix 
T. The antiperiodic boundary conditions are represented by a spin-flip operator P 
that flips the sign of all spins in a given z-slice. 
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The partition function of the periodic system is given by = TrT , while 
the partition function of the antiperiodic system is given by = TrT*P. The 
operators T and P commute and thus have a common set of eigenfunctions. Say 
the eigenvalues of T are Aj and those of P are Pi. The possible values of Pi are 1 
and —1. The partition functions take the form Zp = J2i K = J2i KPi- 

Let us consider the ratio of the partition functions in the broken phase of the 
model. Here the two largest eigenvalues Aqs and Aoa are much larger than the 
other eigenvalues [0 . (The subscripts s and a label eigenfunctions with p = 1 and 
p = —1, respectively.) Hence the ratio of the two partition functions can be well 
approximated by 



(6) 



The corrections are of order {Xis/ Aqs)*- This means that the extension t of the lattice 
in 2;-direction has to be much larger than the bulk correlation length ^. For ,^oa » t 
(Coa = ^Vl'^(i^) is the tunneling correlation length) we can write 

Za t AoaN /„N 

Notice that within this approximation the derivative of G with respect to f3 does 
not depend on t. 

According to this discussion, if ^ « t « ,^oa then already for finite L a surface 
free energy is rather well defined by 

-G + lnt. (8) 

Phenomenologically one can interpret this situation as follows: 
The lattice is short enough that the creation of interfaces in the system with periodic 
boundary conditions is sufficiently suppressed while for the system with antiperiodic 
boundaries only the interface induced by the boundary conditions is present. On the 
other hand the extension of the lattice is large enough not to restrict the fluctuations 
of the interface. 

In order to discuss the L dependence of the surface free energy a model for the 
surface is needed. 

Theoretical studies of the interface are based on the SOS (solid-on-solid) approx- 
imation which essentially neglects overhangs and bulk fluctuations. SOS models 
predict that the roughening transition is of the Kosterlitz-Thouless type [^]. That 



the Ising interface at the roughening point is indeed in the same universality class 
as various SOS models was demonstrated by Hasenbusch using a renormalization 
group matching procedure [jl| . For large /3 there are also rigorous results from linked 



cluster expansions. Borgs and Imbrie have shown [|I^ for sufficiently large /3, i.e., 
when the interface is smooth, that 

Fs^aL\ (9) 
It is believed that this result holds for all (3 > f3r- 
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A model for the interface in the rough phase is the capillary wave model . In 
its quadratic approximation it essentially states that the infrared fluctuations of the 
interface are massless Gaussian. This assumption has been verified numerically in a 
number of cases, see e.g. 0. The massless Gaussian dynamics leads to the following 
finite L behavior of the surface free energy in the rough phase |T^, [TP|, pT| : 



Fs-Cs + aL\ (10) 
Gelfand and Fisher predicted in addition a logarithmic dependence of the surface 



free energy on the lattice size. They did not take into account a prefactor L in the 
partition function that arises when the average position of the interface is fixed via 
a (5-f unction. 

At the roughening transition one has still a Gaussian fixed point, however, with 
logarithmic corrections. Hence eq. (|T^) should be still valid for sufficiently large L. 

It is difficult to compute free energies directly by Monte Carlo (however, cf. [pJ]) ). 
But the derivative of G with respect to /5 is a quantity that can be computed by 
Monte Carlo: 

^ = {H)^-{H),^E.,. (11) 
G can then be obtained by integration over (3: 



G{P) = G{Po)+ f d(3' Es{(3') , (12) 



where /3o is arbitrary. In the case that there is only one interface in the system, 
is the surface energy. The surface energy per area is defined as 

ts = EJL^ . (13) 

In we used the method of 'integration over to determine the surface free 
energies for a wide range of temperatures, for L = 8, 16, 32, 64. We found that the 
surface free energy can be fitted accurately with eq. (plQl). We thus identify the 
coefficient a in front of the factor with the surface tension. 

ASOS Model 

The ASOS model is the solid-on-solid approximation of an interface of a 3D simple 
cubic lattice Ising model (on a (OOl)-lattice plane). It lives on a two-dimensional 
square lattice of size L by L. The partition function is 

Zasos = Z!^^P ' (14) 

h \ <ij> J 

with hi integer. We shall study the quantity 

eAsos = Easos/L^ =< \hi - hj\ > / . (15) 
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In the ASOS limit of the Ising model the inverse temperatures of the two models 
are related by 

Asi„g=/5A505/2. (16) 

In the same limit, there is the following relation between the surface energies: 

Es,lsmg = 2 + 2 Easos (17) 

for corresponding /3 values. The most recent estimate for the roughening coupling 
of the ASOS model is (3r^Asos = 0.8061(3) H]. 



3 Monte Carlo Algorithms 
Ising Model 

Cluster Algorithm for the Ising Interface 

For the production of the Monte Carlo results of ref. and for part of the new 
results to be presented below we employed the Ising interface cluster algorithm of 
Hasenbusch and Meyer p4|. A detailed description of this algorithm can be found 
in 1^. There we computed besides the surface energies also nonlocal quantities like 
the surface thickness and block spin correlation functions. The cluster algorithm 
proved to be very efficient for this purpose. 

Local Demon Algorithm for the Ising Model 

When the focus is on the determination of the energy with antiperiodic boundary 
condition, it is helpful to use a local algorithm instead of the cluster algorithm. It is 
much easier to adapt a local Monte Carlo algorithm for vectorization, parallelization 
or multi-spin coding. 

It turned out that the energy of the system with antiperiodic boundary conditions 
does almost not couple to the slow modes of the local algorithm. 

For the update of the 3-dimensional Ising model we therefore also used a micro- 
canonical demon algorithm I^S], ^ in combination with a particularly efficient 



canonical update pSf of the demons. The algorithm is implemented using the multi- 



spin coding technique 0. Every bit of a computer word carries one Ising spin. 
In order to avoid restrictions of the geometry of the systems 32 independent systems 
are run in parallel. The demon carries three levels carrying the energies 4, 8 and 16. 
The number of demons is chosen to be equal to the number of lattice sites. 

The simulation is done by performing a cycle of 5 groups, where each group 
consists of a microcanonical update of the spin demon system and a translation of 
the demon layer with energy 4,8, or 16 (alternating). Each group is finished by 
updating the demons with energy 4. The hole cycle is completed by updating the 
demons with energy 8. 

The canonical update of a demon layer consists of the following steps 



1. calculate the total number No\d of demons in that layer carrying energy 
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2. calculate a new total number A^new of demons with energy according to the 
probability pe = nEexp{—i3E)/ZD, with = Y,e^eGW{~PE), where ue 
is the number of demon states having the total energy E 

3. if Nncw > Noid then pick randomly A^new — A^oid empty demons and fill them, 
else pick randomly A^^om — A^new occupied demons and clear them 

This is a valid update of the total demon system that provides new total energies 
with a heat bath distribution. Random numbers are only needed for the selection 
of the [A'oid — A'newl demons which change their energy and one for the selection of 
the new total energy. Notice that jA'oid — A^ncwl is only of the order of the square 
root of all the demons. 

Combined with shifts and translations of the demons such an update of the 
demons should be almost as good as a heat bath for every single demon. 

In table |l] we quote the performance of our algorithm on various workstations. 
For comparison we also cite the performance data of algorithms of other authors on 
different supercomputers in table |[ 



VMR Algorithm for the ASOS Model 

The results for the energy of the ASOS model cited in this paper were obtained 
with the help of the VMR (valleys-to-mountains-reflection) algorithm introduced in 
TM. The basic idea of this algorithm is to define valleys and mountains by cutting 



the configuration of height variables with a reflection plane. The flip of a cluster of 
spins can be regarded as reflecting a valley to a mountain (or vice versa). This type 
of algorithm has proved very successful for various SOS models. 



4 Discussion of Results 

Ising Model 

The Surface Energy 

We first explain how we obtained estimates for the surface energy defined in eq. (pT]) . 
The quantity {H)a was always computed using either the interface cluster algorithm 
or the local demon algorithm. For the determination of the energy with periodic 
boundary conditions we used a Fade approximant of the low temperature series for 
the energy (for details cf. 0). The series was put to order 24 in m = exp(— 4/5) by 
Bhanot et al. |2^ and a little later to order 32 by Vohwinkel |3^. In 0], we found 



by comparing with Monte Carlo results, that the use of the Fade approximation was 
safe for (3 > 0.26 for L = 8, /5 > 0.24 for L = 16 and L = 32, and for /3 > 0.235 
for L = 64. For the range of couplings that we want to focus on in this paper 
{P > 0.35), the use of the Fade approximant is safe (for the accuracy required). 
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When computing the surface energy per area one has to deal with systematic 
effects from the finite extension of the lattice in t- and in L-direction.0 

First we carefully studied the t-dependence of the surface energy for three dif- 
ferent /3- values, namely /3 = 0.45 (which is in the smooth phase of the interface), 
/3 = 0.4074 (the roughening point) and f3 = 0.3500 (which is in the rough phase 
of the interface, but still far away from the bulk critical point). The results for 
for interfaces with extension L = 4, 8, 16, 32, 64 and for lattices with various values 
of t are quoted in the tables ^ |^ and ^. The statistics is as follows: For the runs 
with the cluster algorithm we performed 400 000 measurements, separated by an 
update step made up from generating and flipping 8 clusters of two different types, 
cf. 0, and a subsequent Metropolis sweep. For the runs with the demon algorithm 
we made 100 000 measurements on each of the 32 independent copies, separated by 
always 5 cycles as described in section ^ For some sets of parameters we made runs 
with two different random number generators and found perfect consistency. As a 
further check we compared Monte Carlo results for L = t = 4 with the exact result 
for the energy. 

An inspection of the tables reveals that (for fixed L) the estimates for are 
consistent within error bars for t > A and (3 = 0.45. For (3 = 0.4074 we find 
consistency for t > 5 if L < 32 and for t > 6 if L = 64, cf. the discussion of 
finite t effects in section ^ In this context it might be interesting to note that 
the bulk correlation length ^ which governs the exponential decay of the connected 
correlation function is 0.2809 for p = 0.45. For p = 0.4074, one has ^ = 0.3162, and 
for f3 = 0.35 the correlation length is = 0.3897. These estimates are based on a 
Fade evaluation of a low temperature series by Arisue and Fujiwara |31]. Correlation 
length estimates from the series truncated at 13*^* order for (3 = 0.35 differ from the 
Fade approximants in the third digit. 

Using the safe values of t found for the smaller /5- values we computed the interface 
energies presented in table ^ for the range from f3 = 0.35 to (3 = 0.6 in steps of 0.01 
or 0.005. The table also contains estimates for the L = oo interface energy from the 
low temperature series, truncated at 17*^* and at 12*^^ order, cf. the discussion below. 

Concerning the convergence of the surface energies to the infinite L limits we 
make the following observations: For (3 > 0.45 we find that (within the statistical 
accuracy obtained) the estimates converged well, consistent with exponentially small 
L-corrections. Of course, the convergence becomes better for larger (3. 

In the rough phase the interface energy is expected to behave like A + B L^. 
We fitted the /? = 0.35 data for L > 8 and t = 13 with this ansatz and found 
A = -5.96(47) and B = 3.70428(9) with xVdof = 1.135.g 

At the roughening transition the situation is more difficult. We studied the 
differences Es{2L) — Es{L) for L < 64 and found that this quantity scales down 
with a factor of at least 3 always when L is doubled. This gives us confidence that 
Es{oo) - Es{l2%) is not bigger than E,(128) - ^^(64) = 0.0022(1). 

^For a study of finite size effects in a previous high precision simulation performed by Ito see 


^Infinite L-extrapolations according to the law A + B for data in the rough phase were also 



done by Ito 1 12 
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The estimates for the interface energies can be compared with results from the 
low temperature expansion of the same quantity. The series for the interface tension 
is 

17 

a = 2P+Y,anu'' + 0{u^^) . (18) 

n=2 

The coefficients up to order in the low temperature variable u = exp(— 4/3ising) 
were first determined by Weeks et al. They can be found in a paper by Shaw and 
Fisher [|]. The higher coefficients (order through m^^) were computed recently 
by Arisue . The coefficients a„ are quoted in table In order to get from a, 
one has to differentiate eq. ( p^Sf ) with respect to p. The result is a power series in 
u (there is no longer a term ln(M) present). We compared the truncated series and 
its Fade approximants with the Monte Carlo results as follows: For order k, with 
7 < /c < 17, we plotted the result of the truncated series (truncated at order k) 
together with the results of 4 or 5 close-to-diagonal Fade approximants [m/n], with 
m + n = k. For even k we took the five Fades with m = k/2 — 2, . . . , k/2 + 2, for odd 
k we took the four Fades with m = {k — l)/2 — l, . . . ,{k — l)/2 + 2. In all figures the 
truncated series estimates are plotted with a '+' and connected with a broken line. 
For the Fade estimates we use diamonds. On the right hand side of the plots we 
present our Monte Carlo estimates of the surface energy for different lattice sizes. 
The L = 64 result is also plotted with two vertical lines for easier comparison with 
the series results. 

Figure 1 shows this comparison for (3 = 0.50 which is in the smooth phase 
of the interface and far away from the roughening transition. For order > 16 the 
truncated series and the Fade approximants have become consistent with the Monte 
Carlo estimate. Note however, that there is no apparent advantage in using the Fade 
approximants instead of the truncated series. 

f5 = 0.45 is still in the smooth phase. The comparison is summarized in figure 
2. Notice the much larger scale of the y-axis in this plot compared to figure 1. 
The Fade approximants scatter a lot, especially around order 13 where the series 
changes its sign. If one looks only at order < 12, the truncated series seems to be 
even superior to the Fade approximations. 

The scenario becomes even more drastic if one proceeds to the roughening re- 
gion. In figure 3 we show the comparison of the different approximations at the 
roughening point (3 = 0.4074. Here obviously neither the truncated series nor the 
Fade approximants lead to a reasonable approximation. 

The Surface Tension 

The estimates for the surface tension quoted in this paper were obtained with the 
method described in 

However, for the /3-range above 0.35 we used the new and much more precise 
estimates for the surface energy as quoted in table From the results for the 

■^We plan to give a comparison of the different methods to determine the interface tension (see, 
e.g., |, |, |l^, IH) in a forthcoming pubhcation [|| 
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surface energy we determined the surface tension using the method of 'integration 
over /3' as outhned in section |[ In we used two different starting points jSo for 
the integration, namely Po ~ Pc and Pq = 0.6. Both methods yielded compatible 
results. In table ^ we quote our new results obtained by starting the integration 
at the following /? values: For L = 8 and L = 16 we took Pq = 0.545, for L = 32 
we used Po = 0.515, and for L = 64 we started the integration at Po = 0.495. The 
starting point was determined such that the Monte Carlo surface energy estimate 
and the 17*^ order low temperature series for the same quantity are consistent within 
the present error bars. Note that the errors quoted in table |^ are statistical errors 
(lex error bars) that do not include systematic effects. The estimates rely on fits of 
the finite L behavior of the free energy with the law Cs + crL^. 

For the estimation of systematical errors we used the following procedure: One 
defines 

a{L)=F,iL)/L^. (19) 

By the very definition this quantity converges to the interface tension in the infinite 
L limit, however, with stronger finite size effects than definition eq. (plQl). So looking 
at the variation of cr(L) gives one a feeling of the maximal systematic error possible. 

It is also instructive to obtain estimates for a based on the law eq. ([T0|), however, 
using just pairs of adjacent L- values. Then no fit is needed. In table | we quote 
these quantities for P = 0.402359 (which corresponds to m = 0.2) and for P = 0.45. 
We adopt the following rule for the estimation of a systematical error: Take the 
estimates for a from the pair L = 16, 32 and from the pair L = 32, 64 and compute 
the difference. Take this as the systematic error of the interface tension determined 
with the fit method. For the two examples studied in table H we conclude that for 
P = 0.45 the systematic error is smaller than the statistical error. For P = 0.402359 
we arrive at a/{2p) = 0.84487(3). 

In figures 4 and 5 we show the comparison of the Monte Carlo results for the 
two P values quoted above^ with the truncated series and the Fade approximations. 
The Fades were not performed directly for the series for a but (as in |Q) for the 
quantity Q{u) = uexp{2a). 

The conclusions are similar to the ones for the surface energy. Figure 5 demon- 
strates the trap one can get into when the series is too short. Shaw and Fisher 
might have been misled by the convergence and consistency of the Fade and differ- 
ential approximants at order 9 and concluded that the interface tension could well 
be approximated by Fades of ninth order for temperatures below the roughening 
temperature. The now longer series shows that this is a wrong conclusion. Note 
that the Fades seem to converge again at the by now highest available order. But 
still, the value is definitely off from our Monte Carlo estimate.^ 

''We chose u — 0.2 (/3 = 0.402359) for easier comparison with the work of Shaw and Fisher. We 
compared our Monte Carlo results with the series extrapolations also at the more recent estimate 
/3r — 0.4074, with the same conclusions 

'^A preliminary Fade analysis for the interface tension was already performed by Arisue in Q. 
His results are perfectly consistent with ours. However, he could not compare with an independent 
Monte Carlo result 
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Shaw and Fisher pointed out that the use of inhomogeneous differential approxi- 



mants [33| might be superior to using Fades. In fact, these approximants generahze 
Fades and are suitable to deal with functions that have a critical behavior like 
~ A{x){x — Xc)~'^ + B{x). Note, however, that the singularity of the free energy of 
the Ising interface is not of this type. 

Table [1^ shows the results obtained by evaluating inhomogeneous differential 
approximants for a/ {2(3) at /5 = 0.402359 which corresponds to u = 0.2. Like Shaw 
and Fisher, we computed the approximant for the quantity Q{u) defined above and 
then took the logarithm. We here only discuss the order 9 and order 17 approxi- 
mants. The order 9 approximants are fairly stable, however yield too small results 
(as the Fade approximants of this order do). Recall that the Monte Carlo estimate 
for this /5-value is a/{2[3) = 0.84487(3). The order 17 approximants are also fairly 
stable, however, they now overshoot the Monte Carlo estimate definitely. We con- 
clude that using inhomogeneous differential approximants does not cure the problem 
(as was to be expected). The analysis of the surface tension at (3 = 0.4074 leads to 
a similar result. 

In fig. 6 we compare our results for the surface tension with the whole range 
approximant of Fisher et al. ^ in the /5-range 0.35 — 0.50. The approximant 



provided by Fisher and Wen [3l| is obtained as discussed in ref. 0]. In addition the 



critical amplitude of the surface tension is fixed to the value given by Mon |^ . The 
mismatch of the two curves can be explained by the failure of the approximants to 
the low temperature series of order 9 at m = 0.20. The surface tension at u = 0.20 is 
underestimated and since the surface energy is overestimated the gap between the 
two curves increases with decreasing /?. 

For /^-values close to the bulk critical temperature the Monte Carlo result and 
the interpolation are consistent again. This fact confirms the validity of the result 



for the surface tension amplitude obtained by Mon ^6 



The Specific Heat 

The specific heat of models undergoing a KT phase transition was investigated in 
a number of Monte Carlo studies. For the XY model a peak of the specific heat is 
found, which is located in the massive phase of the model at about 0.9/3c, see e.g. 



refs. ^ Swendsen [^] also found a peak of the specific heat for the DGSOS 
and ASOS models in the massive (smooth) phase of the models. For the BCSOS 
model (or F-model) that is an SOS model with the constraint that two neighboring 
height variables must differ by +1 or -1, the specific heat can be computed exactly 



(chapter 8 in the book of Baxter [41 1). The peak of the specific heat lies clearly in 
the massive phase. 

In contrast to these findings, Shaw and Fisher ^ arrive at the conclusion (based 
on their series analysis) that the peak of the specific heat of the 3D Ising interface 
is located in the massless (rough) phase. 

The position of the specific heat peak is, of course, not a universal feature of 
the KT transition since the specific heat stays finite for all temperatures. However, 
the ASOS model is supposed to approximate the Ising interface at the roughening 
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transition even quantitatively quite well. Hence it would be quite surprising if the 
specific heat peak of the Ising interface should be in the massless phase while that 
of the ASOS model is located in the massive phase. 

We computed the derivative of the surface energy with respect to the inverse 
temperature (3 from finite differences of the energy. The results for L = 8, 16, 32 and 
64 are given in fig. 7. For comparison we give the truncated series result for order 12 
and 17. For the lattices of size L = 16, 32 and 64 the derivative of the energy clearly 
exhibits a peak for P ~ 0.43 > Pr- The peak of the specific heat C = itself 
is even slightly deeper in the massless phase. This is in contradiction to the Shaw 
and Fisher result. The peak height still increases considerable for L = 64 compared 
to L = 32. This fact indicates that length scales of order 10 lattice spacings are 
largely involved in the generation of the specific heat peak. On the other hand the 
correlation length at /5 = 0.43 is finite. It should be of order 100. 

For j3 > 0.46 the derivative of the energy obtained from the 12*'^ and 17*^^ order 
truncated series are close together and reproduce the Monte Carlo result within error 
bars. The 12*^ order truncated series contains only coefficients with positive sign 
and is hence increasing monotonically with increasing u. Obviously it cannot predict 
the peak of the specific heat. The situation is slightly different for the 17*^ order 
truncated series. The curve displays a peak at /5 = 0.3747... deep in the massless 
phase, but obviously wrong. 

One expects the low temperature series of the free energy to converge for (3 > (3r- 
Hence in principle one should be able to obtain the specific heat peak accurately 
from the truncated series of sufficiently high order. However, since length scales of 
more than ten are involved one would have to compute diagrams of this extension. 



ASOS Model 

It is instructive to study the same questions in the SOS approximation. The moti- 
vation for this is twofold: First, the comparison of the results demonstrates that the 
problem studied is indeed due to the Kosterlitz-Thouless nature of the roughening 
transition. Furthermore, e.g. in the ASOS model it is much easier to get precise data 
for large L. Table ^ shows our Monte Carlo results for the energy of the ASOS 
model as introduced in section The data were obtained with the help of the VMR 
algorithm. Every single simulation consisted of at least 20 000 updates. 

In table |^ we quote the coefficients of the low temperature expansion of the free 
energy per area for the ASOS model. The expansion variable is u = exp(— 4/3ising), 
with /Rising = /?Asos/2. The series is 

12 

t^Asos = 2/3 + 5^ a„ + C>(n") . (20) 

n=2 

The coefficients a„ with n < 9 are due to Weeks et al. (cited in 0]). We extended 
this series to order 12, using a technique due to Arisue and Fujiwara [0, ||. Our 



result for the ninth order is at odds with the result cited in [Q, where Og = 3185/3 
is quoted. We used our own value in what follows. 
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It is interesting to note that the series (hke the corresponding Ising series) has a 
characteristic change of sign in the coefficients. Prom the series for the free energy 
one easily obtains the series for the energy per area. Figures 8 and 9 compare again 
our Monte Carlo results with the truncated series and with Fade approximants of 
increasing order. The observations are essentially the same as for the Ising model. 
Since we extended the series beyond the order where the change of sign takes place 
we can observe the Fade approximants becoming unstable around this order. 

5 Conclusions 

In this paper, we presented a comparison of Monte Carlo results for interface prop- 
erties with low temperature series. We took the obvious discrepancy between the 
methods as a motivation to improve confidence in the Monte Carlo estimates by 
providing a detailed study of possible systematic errors. 

The failure of the series approximations to improve the series result compared 
with the truncated series as discussed in this article is not completely unexpected. 
However, our detailed study reveals the seriousness of the problem and calls for new 
approaches to deal with essential singularities in series expansions. 
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Figure Captions 



Fig. 1: Ising surface energy at (3 = 0.5000 from Pade, truncated series and Monte 
Carlo. Truncated series estimates are plotted with a '+' and connected with a 
broken line. The close to diagonal Pade estimates are given with diamonds. The 
data with error bars present our Monte Carlo results for the different lattice sizes. 
The L = 64 estimate is also given by two vertical lines. 

Fig. 2: Ising surface energy at /? = 0.4500 from Pade, truncated series and Monte 
Carlo. The symbols are the same as in fig. 1 

Fig. 3: Ising surface energy at /5 = 0.4074 from Pade, truncated series and Monte 
Carlo. The symbols are the same as in fig. 1 

Fig. 4: Ising surface tension at f3 = 0.4500 from Pade, truncated series and Monte 
Carlo. The symbols are the same as in fig. 1 

Fig. 5: Ising surface tension at u = 0.2 [jS = 0.4204) from Pade, truncated series 
and Monte Carlo. For comparison we also quote here the Shaw and Fisher estimate 
(based on the ninth order series evaluation). The symbols are the same as in fig. 1 

Fig. 6: The surface tension obtained by Monte Carlo simulation (dotted line) is 
plotted as a function of the inverse temperature f3. The error bars are smaller than 
the line width. For comparison we give the interpolation result of Fisher et al. 
3^ , |, (solid line). The roughening transition at (3r = 0.4047(3) is indicated by 



a vertical dashed line. 

Fig. 7: The derivative of the surface energy with respect to P is plotted as a function 
of the inverse temperature (3. The dotted lines gives the derivative obtained from 
finite differences of the surface energy for L = 8,16,32, and 64, computed by Monte 
Carlo. The height of the peak grows with increasing L. For comparison we give the 
result from the series truncated at order 12 (Seriesl2) and truncated at order 17 
(Seriesl7). The roughening transition at jSr is indicated by a vertical dashed line. 

Fig. 8: ASOS energy at /5 = 0.85 from Pade, truncated series and MC. The symbols 
are the same as in fig. 1. 

Fig. 9: ASOS energy at f3 = 0.81 from Pade, truncated series and MC. The symbols 
are the same as in fig. 1. 
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machine 


algorithm 


mspins/sec 


SUN ELC 


mc 


2.3 


SUN ELC 


c 


1.1 


SPARC 10 


mc 


7.1 


SPARC 10 


c 


3.2 


HP 9000/735 


mc 


28.0 


HP 9000/735 


c 


14.7 



Table 1: Performance of the multi-spin program on various work- 
stations, given in units of 10^ spin updates per second, 'mc' means 
micro-canonical run, 'c' with canonical update of the demons 



author(s) 


year 


machine 


mspins/sec 


Wansleben et al. |42] 
Creutz et al. [^] 
Kikuchi and Okabe |43| 
Ito and Kanada [^] 

Ito Q 


1984 
1986 
1987 
1990 
1991 


CYBER 205 (two pipes) 
CYBER 205 (two pipes) 
NEC SX-2 
HITAC S820/80 
VP2600/10 


21.2 
117.0 
251.0 
1960.0 
2190.0 



Table 2: Performance of the algorithms of other authors on supercom- 
puters for comparison with table 



t 


L = 4 


L = 8 


L = 16 


L = 32 


L = 64 


2 


2.67433(63) 


2.78942(42) 


2.83759(30) 


2.85006(21) 


2.85225(12) 


3 


2.66940(57) 


2.75544(34) 


2.77860(21) 


2.78114(11) 


2.781092(56) 


4 


2.66799(63) 


2.75238(38) 


2.77445(23) 


2.77717(13) 


2.776936(68) 


5 


2.66775(57) 


2.75246(35) 


2.77440(20) 


2.77677(11) 


2.776790(64) 


9 


2.6700(16) 


2.7517(9) 


2.7746(5) 






10 


2.66645(75) 


2.75223(42) 


2.77442(25) 


2.77683(15) 


2.776835(83) 


17 


2.6688(19) 


2.7541(11) 


2.774(1) 







Table 3: Ising interface energy at (3 = 0.4500 
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t 


L = 4 


L = 8 


L = 16 


L = 32 


L = 64 


2 


2 9542('7l 


3 2408('5l 


3 42063(27) 


3 48668(12") 


3 492507(59") 


3 


2.9479(7) 


3.14030(37) 


3.21718(22) 


3.24322(12) 


3.25120(7) 


4 


2.9419(8) 


3.1244(4) 


3.19578(24) 


3.21901(14) 


3.226227(83) 


5 


2.9414(7) 


3.12342(37) 


3.19428(22) 


3.21745(12) 


3.224334(62) 


6 


2.9427(8) 


3.1230(5) 


3.19410(24) 


3.21735(14) 


3.224389(83) 


9 


2.9384(20) 


3.1235(11) 


3.1936(6) 


3.21723(30) 




10 










3.224219(88) 


11 


2.9520(10) 


3.1237(5) 


3.19443(28) 


3.21756(16) 




17 


2.9431(25) 


3.1247(14) 


3.1941(7) 






33 






3.1946(9) 







Table 4: Ising interface energy at /3 = 0.4074. We made an additional 
run with L = 128 and t — 11. The result for the surface energy is 
3.226375(45) 



t 


L = 4 


L = 8 


L = 16 


L = 32 


L = 64 


2 


3.3173(7) 


3.69719(30) 


3.77344(14) 


3.77574(7) 


3.775673(34) 


3 


3.3684(8) 


3.70942(48) 


3.85405(31) 


3.92648(21) 


3.95577(11) 


4 


3.3479(9) 


3.62530(50) 


3.70230(29) 


3.72218(16) 


3.726973(76) 


5 


3.3432(9) 


3.61092(48) 


3.68357(28) 


3.70146(14) 


3.705558(61) 


6 


3.3437(11) 


3.60982(55) 


3.68145(28) 


3.69886(16) 


3.703239(76) 


7 


3.3446(13) 


3.61025(63) 


3.68106(32) 


3.69870(17) 


3.702884(85) 


9 


3.3431(28) 


3.6094(15) 


3.6816(8) 


3.69854(39) 




13 


3.3440(16) 


3.61122(78) 


3.68080(39) 


3.69851(20) 


3.70282(10) 


17 


3.3411(38) 


3.6102(19) 


3.6836(10) 







Table 5: Ising interface energy at = 0.3500 
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B 


L — H 


L = 16 


L = 32 


L = 64 


)S'i7 


<S'l2 


0.350 


3 61122f78') 


3 68080('39) 




ZU J 


3 70282(11) 


4 59933 

^ • vy ly ty vy vy 


5.77307 




5S'^7Ql'7n'l 


fin45n(''^fi'i 


3 621 57 


'1 9) 


fi9574(Q) 


4 470Q'^ 


5 OQCi"^? 

O .\Jcf\JO i 


0.370 


3 45279f67) 


3 52385f35) 






3 54621(9) 

• ^ vy ^ \ (_y i 


4.23424 


4.56179 


380 


3 3fi878l'fi4'l 






'1 fi) 


4fi4fi9(8) 


Q7'^51 

O . c* 1 OvJ i 


4 14707 


0.390 


3 28038 fSQ") 


3 35522(34) 


O.O 1 oou 


'1 R"! 


3 38049(7) 


3 72531 


3.81746 








(J . (JO -L VJ(J 


'1 6) 




vy .U ±VJUvJ 


fi7785 


0.400 


3 19154('55') 


3 26458('30') 






3 29260(8) 

(_/«^ty^wvyVLy i 


3 50320 

vy • vy vy vy ^ vy 


3 55224 

vy • vy vy ^ ^ 


40^ 








'1 5) 




3 40^00 

vJ.T:V7vJWty 


4'^88Q 

vy .TivJOOty 


0.410 


3 09971 f53'l 


3 16889('29') 


1 Q91 7 


'1 fi\ 


3 19938(9) 


3.31017 


3 33632 

vy . vy vJ VJ vJ^ 






1 1Q79('9Q'l 




'1 fi) 




S 9941 5 


94'^9fi 


0.420 


3 00972f50l 


3 06903(29^ 


OQOOO 


'1 7) 


3 09669(11) 

t_/«vytyvyvyiyv j-J- i 


3.14463 


3 15860 

vy • i vy vy vy vy 


49^ 




3 018nfi('98') 




'18) 




3 071 1 7 


08140 

vy .UO XtiU 


0.430 


2 91929f48) 


2 96759('27~) 

B vy 1 ^_y \ ^ i / 


9 0890"^ 


'1 S'i 


2.98512(11) 


3 00334 

vy • vy vy vy vy 


3 01083 

vy • vy -±- vy vy vy 






9 Q1fifi4('98') 


2 Q2813 


18) 




9 q40fiQ 


9 Q4fi1 8 


0.440 


2 83342 ('47') 


2 86806('27) 

^ • vy v_y \j \_y \ ^ i / 


Z.o ( OUZ 


'1 fi^ 


2 87575(9) 


2 88278 

^•vy^^^ 1 vy 


2 88681 

* vy vy vy vy -l 








9 89450 


'1 fi) 








0.450 


2 75223f42l 


2 77442 f 25) 


2 77682 


15) 


2 77683(8) 

ij • 1 1 vy w vy y Ly 1 


2 77965 

^ • 1 1 ty vy vy 


2 78182 


0.455 




2 73093f23l 




'14) 




2 73370 

^ • 1 vy vy 1 vy 


2 73529 

^* 1 vyvy^ty 


n 4fin 


9 fi7fi88('41 


9 fi8Q51 O'W 




'14) 


9 fiSQQI (7) 


9 fiQI 07 


9 fiQ994 


4fi5 




2 65063('2Tl 




'1 '^) 




2 651 48 


9 65934 


0.470 


2 60753 f38') 


2 61440(20) 


9 fil49fi 


'13) 

1"J7 


2 61408(6) 


2 61465 

^ • Vy -L^Vy vy 


2 61528 

^ . vy X vy Vy 


0.475 






9 ^7Q7Q 


'19) 




2 58036 

ill . VJ OVJ vJVJ 


2 58083 

ill . vJ vy v7 vy vy 


0.480 


2 54599('36') 


9 54823('19) 


9 5481 


'19) 


2 54813(6) 


9 54840 

^ • Vy^VJTiVy 


2.54874 


0.485 








'1 1 ) 




2 51856 

^ • vy -I- vy vy vy 


2 51882 

t vy J- Vy vy ^ 


4Q0 


9 48Qq9('S9') 


9 40041 M 8) 


9 4Qnfin 


'10) 


9 40054(5) 


9 4Q0fi8 


9 4Q087 


5on 


9 44055 fsn"! 


9 44059(1 8) 


9 44014 


'10) 


9 4401 4(5) 


2.44018 


9 44098 


0.510 


2 3961H'28'l 


2 39613(18) 


2 39578 

^ • ly ty 1 


'9) 




2 39580 

^ • vJ iJ v^ O vJ 


2 39586 

^ • vJ tJ v^ vy VJ 


0.520 


2.35725(25) 


2.35658(17) 


2.35657 


;8) 




2.35665 


2.35668 


0.530 


2.32212(23) 


2.32219(16) 


2.32208 


;8) 




2.32198 


2.32200 


0.540 


2.29135(23) 


2.29123(14) 


2.29129 


7) 




2.29118 


2.29119 


0.550 


2.26399(20) 


2.26403(14) 


2.26363 


7) 




2.26374 


2.26375 


0.560 


2.23930(20) 


2.23943(13) 


2.23904 


7) 




2.23922 


2.23923 


0.570 


2.21751(19) 


2.21728(12) 


2.21724 


;6) 




2.21726 


2.21726 


0.580 


2.19748(19) 


2.19757(11) 








2.19755 


2.19755 


0.590 


2.17986(17) 


2.17989(11) 








2.17981 


2.17981 


0.600 


2.16378(17) 


2.16387(10) 








2.16383 


2.16383 



Table 6: Ising interface energy, Monte Carlo and truncated low tem- 
perature series results {Si corresponds to truncation at i*^ order) 
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P "1 K'^P) ^ 


P " 1 K^P) ^ 


.225 .01629(30) 2.3 


.315 .575804(38) 7.0 
^7'^^^94(''^8'l ^ R 


.405 .8499891(27) 6.2 
84QQQ7Q('9fi'l 91 


.230 .04929(25) 4.8 


.320 .597570(31) 6.9 

'^07^8^/' '^^'\ ^ 
.oy / (jO(j\^Oij j u.(j 


.410 .8592883(27) 2.1 

8^Q9Qnfi('99'\ 9 


.235 .08565(18) 1.9 


.325 .618399(28) 5.6 

.U-LO^ZUt Z(J ) O.l 


.415 .8680694(26) 2.0 

8fi8nfi4f!/'99'l 7 8 


.240 .12262(18) 1.1 


.330 .638353(26) 7.4 

.UOOOe/Ul ZO J j.±.o 


.420 .8763354(26) 10.0 
876^93'^('91 4Q 7 


.245 .15971(15) 1.4 


.335 .657448(21) 4.9 

.UO ( (JUUt -LU ) ±0.1 


.425 .8840889(24) 19.2 
884071 4^ Q^l 117 


.250 .19623(14) 1.2 


.340 .675699(13) 2.8 
675751 (^ '?\ 9^ 6 


.430 .8913357(25) 24.3 
8Q1 ^1 54n 8"! 1 78 7 


.255 .23200(12) 1.4 

9'^9n9('19'l 1 


.345 .6931576(53) 3.7 

fiQ'^1 8n9('4Q'l 'iQ f, 
.yjao lOKjz^y^v J ov.v 


.435 .8980930(21) 13.9 

80807*^1 n fi^l IQi^O 


.260 .266799(93) 1.7 


.350 .7098453(66) 1.0 

7(10871 '\((V?'\ '{ I ^^ 


.440 .9043836(18) 5.2 
on r'>riri"iri (V\ 1 7~> 


.265 .300562(97) 1.5 

'^nn^Q9l'Q9'l 9 9 


.355 .7258015(44) 2.5 

79^891 Of'^fi'l 4"^ 


.445 .9102364(18) 2.0 

Q1 0999'^('1 4'\ 1*^7 8 


.270 .333177(90) 2.4 


.360 .7410463(39) 5.6 
741 DfiSI (^A') 111 Q 


.450 .9156833(14) .9 
Q1 f\f^7')f^(^ 9") 07 ^ 


.275 .364607(73) 2.9 

.OU^UOt/lOZ) Z. ( 


.365 .7556038(42) 5.1 


.455 .9207540(12) 1.6 
Q9074fi8('1 9^1 fil 


.280 .394894(71) 1.7 


.370 .7694993(37) 4.7 
7fiQ'i9nfi(''^4'\ fi4 


.460 .9254779(11) 1.6 
.yzu^ 1 ou\^ lu j ou..^ 


.285 .424033(66) 1.8 


.375 .7827565(36) 6.7 
789778n(''^n'\ fi8 ^ 

. J OZ / lOUlOUl uo.o 


.465 .9298820(12) 1.3 
Q9Q878fi('1 0^1 90 


.290 .451998(61) 1.8 


.380 .7953967(37) 8.4 
7Qii41 7A('^'?'\ 7n 4 


.470 .9339903(10) 2.1 

Q^^C188'^('0Q'l Q 9 
.yooyooo\^uyj y.z 


.295 .478878(63) 2.2 

478887('55'l 9 6 


.385 .8074397(33) 9.2 

8074609('3T) 69 4 


.475 .9378262(08) 1.3 

9378951 f07'l 3 9 


.300 .504651(53) 2.8 
.504664(52) 2.9 


.390 .8189052(34) 10.6 
.8189246(26) 66.9 


.480 .9414105(06) 1.0 
.9414100(06) 2.5 


.305 .529376(47) 3.8 

.529393(46) 3.5 


.395 .8298090(32) 10.5 

.8298263(28) 59.4 


.485 .9447624(04) 1.6 

.9447621(03) 1.8 


.310 .553076(44) 5.2 
.553089(42) 4.3 


.400 .8401660(29) 10.1 
.8401801(25) 43.6 


.490 .9478999(04) 1.8 
.9478997(03) 1.1 



Tabic 7: Monte Carlo results for the interface tension. We always 
quote two numbers. For the upper number the L = 16, 32, 64 data were 
included. For lower number also the L — 8 data were included. X 
denotes P^r degree of freedom for the fit of the free energies with 
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P = 0.402359 (u = 0.2) 


L a{L)/{20) 


pair cr/(2/3) 


8 0.849503(17) 
16 0.8460946(87) 
32 0.8451848(42) 
64 0.8449398(30) 


8-16 0.844959 
16-32 0.844882 
32-64 0.844858 


p = 0.45 


L a(L)/(2/3) 


pair a/{2P) 


8 0.916039(10) 
16 0.9156853(56) 
32 0.9156838(24) 
64 0.9156833(14) 


8-16 0.915567 
16-32 0.915683 
32-64 0.915683 



Table 8: Surface tension a{L)/{2(3) together with estimates for the 
same quantity obtained from surface free energies for pairs of L-values 



n 


„lsing 


„ASOS 


2 


- 2 


- 2 


3 


- 2 


-4 


4 


- 10 


- 10 


5 


- 16 


- 24 


6 


- 242/3 


- 194/3 


7 


- 150 


- 172 


8 


-734 


-452 


9 


- 4334/3 


- 3184/3 


10 


- 32122/5 


- 8862/5 


11 


- 10224 


+ 1712 


12 


- 106348/3 


+ 116804/3 


13 


+ 53076 




14 


+ 3491304/7 




15 


+ 74013814/15 




16 


+ 27330236 




17 


+ 160071418 





Table 9: Coefficients of the low temperature series for the 3D Ising 
and ASOS interface tension 
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order 17 


order 15 


order 12 


order 9 


[y/o, 


11 

-LJ 




\7 1^ 11 
['/3> -^1 


n S4Q 1 Q 


[fi/S 11 
[D/C>, iJ 


n 841 Si^ 
U.O^l oo 


[4/9 11 


n 84n 7s 


[O/D, 


11 
J-J 




[D/D, ij 


n R4fi 91 


[^/4 11 


n 84S fi4 


[S/S 11 
[0/0, iJ 


n 84n fiQ 


[7/7 


11 
-■■J 




f^/7 11 




[4/^ 11 


841 9fi 


[9/4 11 




[D/5, 


11 

-■■J 








[S/6 11 


n 841 9^ 






[y/4, 


91 




\7 lA 91 




[^/S 91 




[S/9 91 
[0/Z, ZJ 


n 84n 41 


[0/0, 


91 




[fi/f; 91 

[D/O, ZJ 


n S44 sn 


[4/4 91 




[9/S 91 
[Z/0, ZJ 


n S4n 9Q 




91 




[i^/fi 91 
[D/D, ZJ 


n 84fi 89 


[S/^ 91 
[5/0, ZJ 








[fi/7 


91 


n 84 ii 4^ 


[4/7 91 


n 844 7^ 










r^/8 


91 


n 84^ Q9 
















91 


















"^1 




ffi/4 Si 
[D/4, 




[^/9 Si 
[0/Z, 6\ 


U.O'4:l U ( 


[S/1 Si 
[0/i,OJ 


n S4n 7n 


[D/D, 






[5/5 si 

[0/0, OJ 


n 84 ii 7Q 


[4/S Si 




[9/9 Si 

[Z/Z, OJ 


n 840 0^ 




Si 




[4/6 Si 


845 59 


[S/4 Si 

[0/4, OJ 




[1 /S Si 

[I/O, OJ 


n 840 S9 




Si 








[9/5 Si 

[Z/O, OJ 


849 Ofi 








zLl 




[fi /'^ 41 
[D/ ^, 4J 




[4 /9 41 








[D/0, 


4-1 


RA^ 


[f^/4 41 
[0/4, 4J 




[S/S 41 

[0/0,4] 








f5/6 

\yi 


41 


845 4Q 


[4/5 41 
[^/ ^, ^1 




[2/4 4 


842 08 






[4/7, 


4] 




[3/6,4] 


0.844 90 










[7/3, 


5] 


0.845 63 


[5/3,5] 












[6/4, 


5] 


0.845 68 


[4/4,5] 












[5/5, 


5] 


0.845 63 


[3/5,5] 


0.845 54 










[4/6, 


5] 
















[3/7, 


5] 


0.845 66 















Tabic 10: Surface tension o'/{2j3) at u = 0.2 {(3 = 0.4024) estimated 
from inhomogeneous differential approximants [N/M,L] for Q{u). The 
symbol '-' means that the corresponding approximant could not be 
computed because it did not exist or the numerics was unstable 



I^ASOS 


L = 64 


L = 128 


L = 256 


0.81 


0.70741(57) 


0.70947(33) 


0.70964(18) 


0.82 


0.67948(63) 


0.68052(36) 


0.68106(20) 


0.83 


0.64976(68) 


0.65157(40) 


0.65182(20) 


0.84 


0.64976(68) 


0.62025(37) 


0.62048(20) 


0.85 


0.5874(8) 


0.58837(39) 


0.58869(19) 


0.87 


0.5272(6) 


0.52615(35) 


0.52643(14) 


0.90 


0.4420(6) 


0.44255(26) 


0.44245(14) 


0.95 


0.3320(6) 


0.33177(28) 


0.33173(15) 


1.00 


0.2517(4) 


0.25153(23) 


0.25128(10) 


1.10 


0.1496(3) 


0.14869(15) 


0.14905(7) 



Table 11: Monte Carlo results for the energy of the ASOS model 



21 



2.4415 

2.441 - 

2.4405 - 

2.44 - 

2.4395 - 

2.439 - 

2.4385 - 

2.438 - 

2.4375 - 



o 




Fig. 1: Ising surface energy at beta=0.5000 from Pade, truncated series andMC 
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Fig. 2: Ising surface energy at beta=0.4500 from Pade, truncated series andMC 
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Fig. 3: Ising surface energy at beta=0.4074 from Pade, truncated series andMC 
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Fig. 4: Ising surface tension at beta=0.4500 from Pade, truncated series andMC 
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Fig. 5: Ising surface tension at beta=0.4024 (u=0.2) from Pade, truncatedseries and MC 
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Fig. 8: ASOS energy at beta=0.85 from Pade, truncated series and MC 
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Fig. 9: ASOS energy at beta=0.81 from Pade, truncated series and MC 
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